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Electronic structure calculations are presented for the well-known CuGe03 and the recently dis- 
covered isostructural CuSiOs compounds. The magnitude of the dispersion in chain direction is 
considerably smaller for CuSiC>3, whereas the main interchain couplings are rather similar in both 
compounds. Starting from extended one-band tight-binding models fitted to the bandstructures, the 
exchange integrals were estimated for both compounds in terms of a spatially anisotropic Heisenberg 
model. Remarkable frustrating second neighbor couplings are found both for intra- and inter-chain 
interactions. A magnetic moment of about 0.35 [ib is predicted for CuSiOs in the Neel state. 

PACS numbers: 71.15-Mb 71.20.-b 75.30.Et 

Low-dimensional spin systems such as chains or lad- 
ders are of fundamental interest for contemporary solid 
state physics due to their peculiar electronic and mag- 
netic properties. During the last years, many related ma- 
terials have been found within the cuprate family, famous 
for the high temperature superconductivity. All cuprates 
contain Cu04 plaquettes. In most cases it is energetically 
favorable to connect these plaquettes by the formation of 
chains or planes. According to the number (n = 1,2) of 
oxygen atoms shared by adjacent plaquettes, these com- 
pounds can be classified as so-called edge-shared (n = 2) 
or corner-shared (n = 1) compounds. 

Obviously, the type of sharing affects strongly the 
physical properties of the compounds under considera- 
tion. For example, corner sharing leads to strong anti- 
ferromagnetic coupling between neighboring plaquettes 
compared with the weak inter-chain interactions.cl As 
a result, the straight CUO3 chain in S^CuOs is the 
best known realization of the one-dimensional spin-1/2 
Heisenberg modelp with an in-chain exchange coupling 
of about 2200 K, but with a Neel temperature of only 5 K 
and with an extremely small ordered magnetic moment 
of about 0.06 //b P both due to a small residual interchain 
exchange coupling. Spin-charge separation in the excita- 
tion spectra could be observed for S^CuOs and for the 
double chain compound SrCuOa-P 

Somewhat surprisingly, in contrast to the similarity 
between different corner-shared chain compounds, the 
magnetic properties in the edge-shared chain family ex- 
hibit a remarkable variance. Thus, the edge-shared 
CuC>2 plaquettes in Li 2 Cu02 order antiferromagnetically 
with a ferromagnetic arrangement along those chains and 
with a large ordered moment of 0.9 whereas the 

same chain in CuGeOa shows a spin-Peierls transition at 
low temperatures.El Antiferromagnetically ordered chains 
were observed in Cug-^Zn^GeOa for small concentra- 
tions of Zn impuritiesB It is noteworthy that, even for the 
intensively studied CuGe03, a consensus with respect to 
the quantitative description of competing or complemen- 
tary interactions such as the inter-chain coupling, frus- 
tration and spin-phonon coupling has not been reached so 




FIG. 1. The orthorhombic unit cell of the CuSiOs-crystal, 
perspective view (top), front view (down left) and top view 
(down right). The edge-shared cuprate-chains run along the 
c direction and are canted against each other. 

far, Hi despite the achieved qualitative understanding of 
their influence on different magnetically ordered states. 

Naturally, the magnetic properties depend very sen- 
sitively on the electronic interactions in these systems. 
Therefore, a comparative study of the electronic proper- 
ties of closely related systems can shed light on the inter- 
actions responsible for the magnetically ordered states 
mentioned above. In this context, the recent discovery 
and firstaavestigations of the long searched for compound 
CuSi03,EJ which is isostructural to the prototypical in- 
organic spin-Peierls system CuGe03 is of great scientific 
interest. The crystal structure of CuSiOs is shown in 
Fig. |]. The most important feature for the magnetic 
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properties are the planar edge-shared CuC>2 chains run- 
ning along c-direction. These chains are very similar to 
those of CuGe0 3 . The Cu-0(2) bond length in CuSi0 3 
(CuGe0 3 ) is 1.941 A (1.942 A), the Cu-0(2)-Cu bonding 
angle is 94° (99°). 

Thus, the question arises, whether the very recently 
observed phase transitiorO near 8 K does point to a 
new inorganic spin-Peierls system or to another ordered 
state realized at low temperature. To get theoretical in- 
sight into possible scenarios, we present here comparative 
band-structure calculations and tight-binding examina- 
tions for CuSi03 and CuGeC>3. In this context we note 
that for the latter compound several (non full-potential) 
bandstructure calculation have been reported (e.g. in 
Rcf. |l2|), but to our knowledge the inter-chain interac- 
tion has not been analyzed in detail. 

The relevant electronic structure of these materi- 
als is very sensitive to details of hybridization and 
charge balance. In order to obtain a realistic and 
reliable hopping part of a tight binding Hamiltonian, 
band-structure calculations were performed using the 
full-potential nonorthogonal local-orbital minimum-basis 
schemeE3 within the local density approximation (LDA) . 
In the scalar relativistic calculations we used the ept-, 
change and correlation potential of Perdew and Zunger.til 
Cu(4s, ip, 3d), 0(2s, 2p, 3d), Ge(3d, 4s, 4p, id) and 
Si(2p, 3s, 3p, 3d) states, respectively, were chosen as 
minimum basis set. All lower lying states were treated 
as core states. The inclusion of Ge 3d and Si 2p states 
in the valence states was necessary to account for non- 
negligible core-core overlaps. The O and Si 3d as well as 
the Ge Ad states were taken into account to increase the 
completeness of the basis set. The spatial extension pf 
the basis orbitals, controlled by a confining potentially 
(r/ro) 4 , was optimized to minimize the total energjt 

The results of the paramagnetic calculation!!!?! for 
CuSi0 3 (see Fig. | (a)) and CuGe0 3 (see Fig. | (b); 
we find similar results as the non full-potential calcula- 
tion of Ref. |l2]) show a valence band complex of about 
10 eV width with two bands crossing the Fermi level in 
both cases. These two bands are well separated from 
the rest of the valence band complex and show mainly 
Cu 3d and 0(2) 2p character in the analysis of the cor- 
responding partial densities of states (not shown). We 
note that the occupancy of the two 0(2) 2p orbitals along 
and perpendicular to the chain (lying in the plaquette- 
planes) is rather different, but it is almost identical for 
the corresponding orbitals in both compounds. Therein, 
we found only a small admixture of 0(1) 2p and Ge 4s 
and Ap states, respectively, with a total amount of few 
percent. The examination of the eigenstates of the lat- 
ter bands at high symmetry points yields an antibonding 
character typical for cuprates. Here these relatively nar- 
row antibonding bands are half-filled. Therefore, strong 
correlation effects can be expected which explain the ex- 
perimentally observed insulating groundstate. Despite 
almost perfect qualitative one to one correspondence of 
all valence bands and main peak structures in the den- 



sities of states (DOS) (compare right panels in Fig. ||), 
the most important differences between both compounds 
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FIG. 2. Band structure and total density of states for 
CuSiOa (a), CuGeC>3 (b), and the zoomed antibonding bands 
(c) (CuSiOs full lines, CuGeC>3 dashed lines). The Fermi level 
is at zero energy. The notation of the symmetry points is as 
follows: Y = (010), T = (Oil), Z = (001), X = (100), S = 
(110), A = (111). The chain direction corresponds to Y-T, 
Z-F and S-A. 

occur for the antibonding bands (shown in detail in 
Fig. H(c)). Therefore, we restrict ourselves to the ex- 
tended tight-binding analysis and the discussion of these 
antibonding bands. 

The dispersion of these bands has been analyzed 
in terms of nearest neighbor transfer (NN), next 
nearest neighbor transfer (NNN) and higher neigh- 
bor terms in chain direction, but only NN hopping 
and a diagonal transition term between the Cu02- 
chains have been considered (see Fig. ||). Then, 
the corresponding dispersion relation takes the form 
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FIG. 3. Schematical chain and stack arrangement of 
Cu02-plaquettes, respectively, and considered transfer pro- 
cesses within the bc-plane (left panel) and in the ab-plane 
(right panel). 



E(k) 



m=lA 



tmz cos(mz) + cos(;r) [t x + 2t xz cos(z)] 



+ cos(y/2) [t y + 2t yz cos(z) + 2t xy cos(x)} ) , (1) 

where x = k z a, y = k y b, z = k z c. Notice that in our effec- 
tive one-band description the upper band (see Fig. ||(c)) 
e.g. along T-X corresponds to k y = 0, whereas the lower 
one corresponds to k y — 2n/b. The assignment of the pa- 
rameters has been achieved by two numerically indepen- 
dent procedures: By straightforward least square fitting 
of the whole antibonding band in all directions and by 
using the bandwidths, the slopes and the curvatures at 
special selected high symmetry points. The latter proce- 
dure has the advantage to be less affected by hybridiza- 
tion effects from lower lying bands near the bottom of 
the antibonding band (being of some relevance near the 
Z-point in Fig. |2j). 
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CuGeOs 


-175 
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-5.5 
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-20.6 


CuSiOs 


-88 


-31 


-4.5 


-2.4 


-36 


-21.2 



TABLE I. Transfer integrals ti (in meV) of the extended 
one-band tight-binding model for CuGeOs and CuSiOs. The 
remaining omitted terms from Eq. (1) were found to be irrel- 
evant. 

The results are shown in Tab. ||. The errors can be esti- 
mated between 1% for the large and 10% for the small pa- 
rameters from the difference of both mentioned above fit- 
ting procedures. The analyzed antibonding bands of both 
compounds exhibit a rather similar shape except near the 
Z-points, where the hybridization with lower lying bands 
produces an additional band-crossing for CuGeC>3 (see 
Fig. ^(c)). Recall that the main difference to the corner- 
shared chains as e.g. in Sr2Cu03 is a much smaller in- 
chain NN transfer due to the different geometry. 

In spite of the qualitative similarity, the calculated val- 
ues for the transfer integrals are quite different. The in- 
chain dispersion is nearly twice as large for CuGeOs in 



comparison to CuSiOs. This can be attributed mainly 
to the larger Cu-O-Cu bond angle in CuGe03 (99° and 
94°, respectively). However, this geometrical effect is 
somewhat reduced by the different on-site energies of the 
oxygen orbitals along and perpendicular to the chain (ly- 
ing in the plaquettes planes). The latter difference is 
reflected by the larger separation of the corresponding 
bands at the Z-point in CuSiC>3 (see Fig. ||). 

The inter-chain dispersions in b direction are compara- 
ble. For both compounds, we find also rather significant 
diagonal hopping terms t yz which are reflected by differ- 
ent dispersions along the X-S and the T-Z directions. 
Somewhat surprisingly, we found a sizeable dispersion in 
x-direction for CuGe03 but only a very weak one for the 
CuSiOs counterpart. 

From the transfer integrals discussed above, we con- 
clude that both compounds are not so well-defined quasi 
one-dimensional systems as compared to the corner- 
shared CUO3 chain compoundsEEZl. The inter-chain cou- 
pling is rather significant for CuGeOs, and CuSiOs can 
even be regarded as an anisotropic two-dimensional sys- 
tem. Since increasing inter-chain coupling tends to desta- 
bilize the spin-Peierls state&3, a Neel ordered antiferro- 
magnetic ground state might be expected for CuSiOs in 
contrast to the spin-Peierls state realized in CuGeOs. 

The obtained transfer integrals enables us to estimate 
the relevant exchange integrals J. This knowledge is cru- 
cial for the derivation and examination of magnetic model 
Hamiltonians of the spin- 1/2 Heisenberg type frequently 
used in the literature: 



ijSiSj 



(2) 



In general, the total exchange J can be divided into 
an antiferromagnetic and a ferromagnetic contribution 
J = J AFM -\- J FM . In the strongly correlated limit, 
valid for typical cuprates, the former can be calculated in 
terms of the one-band extended Hubbard model j{^ FM 
= Atfj/(U — Vij). The indices i and j correspond to near- 
est and next nearest neighbors, U is the on-site Coulomb 
repulsion and Vij is the-inter-site Coulomb interaction. 
From experimental dataMI mapped from the standard pd- 
model onto the one-band description, one estimates U—V 
~ 4.2 eV. For the sake of simplicity, we neglect the dif- 
ference in the quantity U — V in the compounds. The 
calculated values for the exchange integrals are given in 
Tab. 0. 
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CuGe0 3 
CuSiOs 


29 
7.4 


15 
3.8 


2.5 
0.9 


0.4 
0.006 


1.11 

1.25 


0.4 
0.43 


0.17 
0.35 


0.21 
unknown 



TABLE II. Exchange parameters Ji (in meV) for CuGeOs 
and CuSiOs, and local magnetic moments (in /is) in the Neel 
state derived from them (see text). The experimental value 
fi ex is an average over various studies mentioned in the text. 
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The value of the NN exchange integral J^ FM ~ 30 meV 
in CuGeOa exceeds the experimental values r©f about 
11 meV from inelastic neutron scattering dataEa, about 
14 meV from magnetic suseeptibilitytil and about 22 
meV from Raman scatteringca. This points to a signifi- 
cant ferromagnetic contribution due to, the Goodenough- 
Kanamori- Anderson-type interactionEl In the following, 
we shall adopt 15 meV for the resulting total exchange 
coupling Ji as a representative value, suggested by the 
average of the above mentioned experimental data. Ow- 
ing to the lack of experimental data we assume the same 
ratio J\/Ji FM in CuSiO-3 as in CuGe03, suggested by 
the quite similar 0(2) 2p orbital occupancies mentioned 
above. For the latter compound, we note the reason- 
able agreement with the available experimental data and 
most of our calculated antiferromagnetic values for the 
remaining exchange parameters. Hence, further possible 
ferromagnetic contributions seem to be less relevant and 
are neglected in the following considerations. 

Further simplification can be obtained mapping Ji and 
the frustrated NN term J«-.onto an effective intra-chain 
coupling J|| — J 1 — 1.12 J2.BJ The calculated values for Jn 
are 12.2 meV for CuGe0 3 and 2.8 meV for CuSi0 3 , re- 
spectively. The latter value is close to the value of 2 meV 
reported by Baenitz et al. from a one-dimensional fit of 
magnetic susceptibility dataO We find also a consider- 
able inter-chain frustration J yz = j3J y with /3=0.36 (0.34) 
for the Ge- (Si-) compound. [This is in good agreement 
with the suggestions of Uhrigp (3 w 0.5 for CuGe03. 

Transfering the above mentioned idea if^ map frus- 



trating terms onto one effective coupling,cil we adopt 
J± = J y — 2 J yz for the effective inter-chain exchange pa- 
rameters in 6-direction. The factor of two is introduced to 
account approximately the twice as large number of sec- 
ond neighbors. The effective anisotropy ratio R = J±/J\\ 
measures approximately the magnitude of quantum fluc- 
tuations. In the crossover region between one and two 
dimensions, quantum fluctuations do strongly affect the 
magnitude of the staggered magnetization m and the lo- 
cal Cu moment y, = g^udm 21 T = for a Neel ground 
state, where (?l=2.06 to 2.26c3 denotes the (anisotropic) 
Lande-factor (tensor) for Cu 2+ in CuGeOa and rid ~ 0.8 
is the hole occupation number of the related Cu 3d pla- 
quette orbital. Using the expression 



m = 0.39\/fl(l + 0.095i?) ln 1/3 (1.3/i?), 



(3) 



taken frorp_Ref. gfj, we arrive at 0.17/xb in_ceasonable 
agreementcll with the neutron data 0.22±0.02c3 and 0.2cB 
for the disorder induced Neel state achieved below 4.5K 
in Zn-doped CuGc03. The same approach predicts a 
significantly larger value of about 0.35/is for CuSi03 re- 
alized in a possible Neel state. 

To summarize, our LDA-FPLO calculation reveals 
valuable insight into the relevant couplings of CuGe03 
and CuSiOs. We can classify CuGe03 as a quasi one- 
dimensional compound with significant inter-chain inter- 
action, whereas CuSiOs is closer to an anisotropic two- 
dimensional compound. The significantly reduced energy 



scale of the in-chain exchange interactions and the large 
inter-chain interaction in CuSiOs are less favorable for a 
spin-Peierls state than for a Neel order. However, due to 
the large frustrations other states such as a spin-Peierls 
state cannot be excluded. Further investigations are re- 
quired to elucidate the unknown ground state. 

We acknowledge fruitful discussions with M. Baenitz, 
C. Geibel, W. Pickett and G. Uhrig. This work was sup- 
ported by individual grants of the DAAD (H.R.) and the 
DFG (S.D.). 
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